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Abstract. 

Mass calculations carried out by Strutinsky's shell correction method are based 
on the notion of smooth single particle level density. The smoothing procedure is 
always performed using curvature correction. In the presence of curvature correction 
a smooth function remains unchanged if smoothing is applied. Two new curvature 
correction methods are introduced. The performance of the standard and new methods 
are investigated using harmonic oscillator and realistic potentials. 
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1. Introduction 

Masses, or, equivalently ground state energies of nuclei are very important basic 
quantities. For example, shell structure changes around drip lines or existence of new 
magic numbers beyond the heaviest known nuclei can be studied with the help of mass 
measurements. Important applications are the determination of possible decay modes 
of a system and the calculation of energy release in a given reaction. Mass values are 
important not only for nuclear physics but for astrophysics and particle physics as well. 
Nowadays not only the precision of mass measurements has improved dramatically but 
it is possible to measure masses of short lived nuclei [H |2] . Very accurate experimental 
mass values are now available for almost every part of the periodic table and this is a 
great challenge for theoretical description. 

Global mass calculations have been carried out using the so called macroscopic- 
microscopic method [3l H] . In this approach the binding energy is the sum of two terms 
E{N, Z) = Emac{N, Z) + Emic{N, Z) . The first term Emac{N, Z) is a smooth function 
of the neutron (A^) and proton number (Z). It can be determined using a macroscopic 
model. The simplest macroscopic expression is derived from the liquid drop model. This 
form can be generalized in several respects. E.g. a well known macroscopic term is the 
formula of the finite range droplet model 0] . 

Beside the smooth trend the binding energy formula contains an oscillatory part 
with pure quantum origin. The microscopic part Emic{N, Z) consists of two terms shell 
correction and pairing energy Emic{N, Z) = Esheii{N, Z)+Epair{N, Z). The pairing term 
is usually calculated using BCS or Lipkin-Nogami scheme and the shell correction by 
the Strutinsky's method [5]. 

The first step in the calculation of shell correction is to take a single particle 
Hamiltonian. The energy E of the ground state is simply sum of the energies of occupied 
single particle states. This energy can be decomposed into a sum of a smooth E and a 
fluctuating term E, E = E + E. It can be argued using the Strutinsky's energy theorem 
^ that the quantities EgheiiiN, Z) and E are identical. 

Recently new facts have emerged concerning the fluctuating part of the energy. 
The shell correction is a rapidly changing function of the particle number. Earlier, 
the average value of E was expected to be zero. Recently it has been shown that the 
average value of the fluctuating part of the energy is non zero and it contributes to the 
smooth part of the energy [7j. The extra smooth trend of the average value of E can be 
analytically calculated [TJ [8] in case of simple model potentials. It would be desirable 
to minimize this extra smooth trend of the shell correction. A study of the nuclear 
mass error distribution [9] , using results of the Strutinsky type macroscopic-microscopic 
calculation of has been revealed a regular correlation and it has been suggested that 
refinement of the Strutinsky's method is needed. 

The focus of the present paper is how to calculate the smooth part of the energy E 
or equivalently the shell correction. The basic concepts are the quantum single particle 
level density (SPLD) and the smooth SPLD. The quantum SPLD is simply sum of 
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Dirac-delta functions centered at the single particle energies. The smooth SPLD is 
free from singularities and it is a smooth function. It may be determined using semi- 
classical methods ^0\. Shell correction calculations [SI [H] were carried out using the 
smooth semi-classical SPLD based on Wigner-Kirkwood expansion. Smooth SPLD can 
be derived using only the single particle energies and a special smoothing method called 
Strutinsky's smoothing or standard averaging method [11] . The Strutinsky type |1] 
and semi-classical [3] macroscopic-microscopic models differ on the way how the smooth 
SPLD is calculated. Beside the energy averaging particle number space averagings were 
developed [121 [131 [E] in order to calculate E and even the combination of averaging 
spaces was considered [T5] . 

The Strutinsky's smoothing of quantum SPLD is carried out with a well defined 
procedure. First, a normalized, positive and even function is selected. The smoothing 
means that the convolution of the quantum SPLD and the kernel function is calculated 
[6l [16]. This simple procedure however is insufficient. We have to demand self- 
consistency of the smoothing process. If the smoothing is applied to a sufficiently smooth 
function (polynomial) then the result has to be identical with the original function. This 
property is referred to as curvature corrected smoothing. 

In the most standard form of the Strutinsky's method the curvature correction is 
taken into account by multiplying the starting kernel function with an appropriately 
chosen polynomial [6]. Brack and Pauli have introduced [I6] an alternative method 
to achieve the effect of the curvature correction. In this paper we introduce two new 
curvature correction methods. 

The organization of the paper is the following. In section 2 we describe the 
Strutinsky's shell correction method and than in section 3 we give precise meaning 
of the curvature correction. Two standard methods and our new methods for curvature 
correction are also presented here. Finally at section 4 we compare the methods using 
harmonic oscillator and Saxon- Woods potentials. Summary is given at section 5. 

2. Strutinsky's shell correction method 

In the Strutinsky's method first an appropriate single particle Hamiltonian is taken. 
The eigenvalues of bound states and the corresponding degeneracies are denoted by ej 
and di, respectively. For simplicity the {ej} series is arranged in increasing order. An 
important notion of the method is the single particle level density (SPLD). The so called 
exact or quantum mechanical SPLD is of the form 



where 6{e — e^) is the Dirac-delta function and the number of bound states is Nb. 




(1) 



i=l 
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The central notion of the Strutinsky's method is the smooth SPLD. It is denoted 
by g{€). The smoothed Fermi-level A is obtained from the particle number equation 

Nr = f g{e)de, (2) 

J — oo 

where is the neutron or proton number of the considered nucleus since the shell 
correction is separately calculated for neutrons and protons. The smoothed total energy 
E can be calculated by an integral 

E = [ eg{e)de. (3) 

J — oo 

The shell correction is defined by the equation 

E = Esp- E, (4) 

where the total single particle energy Egp is taken in the spirit of the independent particle 
shell model 

i=l 

The new energy series {e-} is derived from the series {cj} by repeating the eigenvalues 
ej di times. 

2.1. Smoothing with convolution 

Convolution is an ideal tool to smooth an oscillatory function y{x). We take a bell 
shaped kernel function K{x) and the smoothed y{x) function is calculated by the 
convolution of K{x) and y{x) 

y{x) = - r K(^—^)y{x')dx\ (6) 

where the positive real number 7 defines the smoothing width. The smooth SPLD can 
be obtained from (E]) by identifying y with the quantum SPLD g. The result is 

-g{e) = -J2d.K('-^). (7) 

The particle number equation takes the form 

Nr = ^di I K{x)dx (8) 
i=i 

and the smoothed total energy is 

E = ^ c/j I 7 / xK{x)dx + ei K{x)dx J , (9) 

\ J ~oo J —00 J 

where i/j = (A - 6^)77. 
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3. Curvature correction 

We expect that a proper definition of smootliing lias tlie following property: the 
smoothing does not change a smooth function. Here we identify smooth functions with 
polynomials [HI E]. We require that the mapping y{x) — )■ y{x) is the identity map on 
the set of polynomials of order at most (2m + 1) i.e. 

Pm{x)=Pm{x) = - K[ )pm{x)dx, (10) 



where Pm{x) is an arbitrary polynomial of order at most (2m + 1). We will say that a 
kernel K satisfying (fTOl) has curvature correction of order m. If we want to emphasize 
this property the kernel will be denoted by instead of K. A kernel with curvature 
correction of order zero will be called starting kernel. If we want to stress actual values 
of the smoothing range 7 and the order of the curvature correction m we will use the 
notation ^(e) for the smooth SPLD. 

3.1. Standard curvature correction methods 

For completeness we review two standard curvature correction methods. In the most 
widely used method a kernel with curvature correction of order m is searched in the 
form [6] 

Kmit) = Pm{t)Ko{t), (11) 

where Pm{t) is a polynomial of order 2m 

m 

P^{t) = Y,Cit'\ (12) 

The coefficients q are determined by the condition ( fTOl) . The following set of equations 
can be gained [6] 

m 

Y,M2ii+j){Ko)c2i = Sj,o 0<3<m, (13) 

i=0 

where the moments of an arbitrary kernel K is defined by 

POO 

Mi{K)= fK{t)dt, i = 0,1,2,.... (14) 



This approach is called polynomial curvature correction (PCC). 

Brack and Pauli [I6] have taken the kernel with curvature correction of order m in 
the form 

i^r.(t) = 5^a2.— i^o(t). (15) 

i=0 

We will refer to this method as Brack and Pauli curvature correction (BPCC) method. 
The condition ( ITOj) uniquely determines the coefficients a2i- The value of oq is 1. The 
other expansion coefficients satisfy the following set of equations 

m 

Bija2j = -h2i (z = 1, . . . m), (16) 
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where Bj 



moments of Kq 
h2i 



1, = (i < ]] 



i2^)\ 



^2{i-^j) (1 < « < i) and is related to the 

(17) 



3.2. Curvature correction with many widths 
The error term of the smoothing 
5y[x) = y{x) - y{x) 
can be cast into the form 



Sy{x) 



oo 
n=l 



n! 



(19) 



where the n-th derivative of y{x)is denoted by y^"\x) and it is assumed Mq(K) = 1. 

If the moments of the kernel function K satisfy the following set of equations 
Mo{K) = 1, Mi{K) = 0, i = 1, 2, . . . , (2m + 1) then from ([19]) it is clear that such a 
kernel function has curvature correction of order m. 

Let us assume that we have a kernel Km and suppose that it is an even function. 
We claim that a kernel with curvature correction of order (m + 1) can be written in the 
form 

Km+l{t) = — 



m > 0. 



Here c is an arbitrary real number but c 7^ and c 7^ 1. 

In order to prove our claim it is enough to show the validity of the equations 



(20) 



and 



Km+iit)dt = l 



fKm+i{t)dt = 0, 



1,2, 



(2m + 3). 



(21) 



(22) 



The proofs are very simple. We have to substitute fl2U]) into fl?Il) and (1221) . use the 
linearity of the integration and in the second term of the left hand sides change the 
integration variable from t to u = t/c. 

If the same c value is used at each iteration step (starting from m = 0) then the 
following closed expression for the kernel Km (m > 1) can be derived 



Km{t) 

where the notation 



A 



m 
k=0 



m 













(23) 



m 
r 



1 

(l-g'")(l-g'"~^)...(l-g'"~''+^) 





if r = 

if 1 < r < m 

if r > m 



(24) 
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stands for the Gaussian coefficient (q-binomial coefficient) and 

^ fc=0 ^ ^fc=l ^ 

Kernels derived using (!20|) or ( 123|) are called many width kernels and the smoothing 
method is named as many width curvature correction (MWCC). The name comes from 
the form of the smooth SPLD of the MWCC method 

m 

g,,M = A^^(-l)V"-'=)("'-'=+2> 

k=0 

The smooth SPLD corresponding to a many width kernel with curvature correction of 
order m is a linear combination of smooth SPLD's derived from the starting kernel 
having smoothing width 70^^ {k = 0, ...m). The individual terms g-yckfii^) behave 
differently, some terms may show the shell structure, other terms may over smooth but 
the proper linear combination can make a balance. 



(26) 



3.3. Derivative curvature correction 

The recursion fl20l) is not valid at c = 1. If the limit c — 1 is taken using the L 'Hospital 
rule then a new kernel with curvature correction of order m emerges 

2m + 3 Id 

Km+lit) = K.Jt) + t — Krnit). (27) 

2m + 2 ^ ^ 2m + 2 ^ ^ ^ ^ 

Kernels obtained in this way are called derivative kernels and the method is called 
derivative curvature correction (DCC). 

With the aid of the recursion fl27|) a closed expression can be derived for the 
derivative kernels (m > 0) 

Km{t) = y ^^^^^4^^t^^i^o(t), (28) 
^ ' ^ m!2'" dt^ ^ ' ^ ' 

k=0 

where the coefficients a{n, m) obey the recursion 

a(n, m) = (2n + m — 2)a{n — 1, m) + a{n — 1, m — 1) (29) 

with initial conditions a(l, 1) = 1, a(n, 0) = (n is arbitrary) and a{n,m) = for 
n < m. Apart from the multiplicative factor t'^ the kernel ( 128|) is similar to the kernel 
( !T5|) . In contrast to the BPCC method the linear combination coefficients in (128|) are 
independent from the starting kernel. In the BPCC method the linear combination 
coefficients are calculated using f lT6|) . 

The recursion ( 127|) simplifies the calculation of auxiliary quantities of the shell 
correction method. For the calculation of the particle number we need the integral (see 

(ED) 

Nm{x) = [ K^{t)dt, m > 0. (30) 
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We can easily get the following recursion for the function Nm{x) 

Nm+i{x) = N^x) + -^—xKmix). (31) 

2m + 2 

This recursion is valid provided limt_j._oo (^Kjn{t)^ =0, m > 0. For the calculation of 
the smoothed energy we need the integral (see Qj) 

E^{x) = I tKm{t)dt, m > 0. (32) 

J — oo 

We can derive the following recursion 

2m + 1 1 

Em+i(x) = Em(x) H x'^Km(x). (33) 

+^ ^ 2m + 2 ^ ^ 2m + 2 ^ ^ ^ ^ 

Here we assumed limt_^_oo (^"^Kmit)^ =0, m > 0. 



4. Results 



At the beginning of an application a starting kernel has to be chosen. The most 
widely used form corresponds to a Gaussian Kf){t) = exp(— t^)/i/7r. The kernel 
Ko{t) = 1/(2 Cosh2(t)) was considered in [16] using the BPCC method. We will 
use both functions and call them Gaussian and Cosh-type kernels respectively. An 
interesting observation have to be mentioned. In the case of Gaussian starting kernel 
the methods PCC, BPCC and DCC lead to the same kernel Km- The MWCC method 
however furnishes new kernel function. In the case of Cosh-type starting kernel each 
presented method is different. 

The MWCC method in addition to the smoothing width 7 and the order of the 
curvature correction m contains a new dimensionless technical parameter c. At each 
fixed values of 7 and m the value of c is determined with the following procedure. An 
averaging region around the investigated particle number A^^- is selected. In the following 
calculations the size of the averaging window was 100. Shell corrections of all nucleus 
belonging to the averaging window are calculated. The c parameter is used to minimize 
the average value of the shell corrections. This procedure is inspired by recent works 
[TJ [8]. The shell correction is a rapidly changing function of the particle number. The 
average value of E was expected to be zero but it turned out that the average of the 
fluctuating part of the energy is non zero and it contributes to the smooth part of the 
energy [7]. The value of c is determined in such a way to minimize this extra smooth 
trend of the shell correction. 

The shell correction method contains two technical parameters the smoothing width 
7 and the order of the curvature correction m. The result is supposed to be independent 
from the values of these parameters if they are selected in a reasonable way. In the 
case of harmonic oscillator or similar potentials (infinite number of bound states) it 
was demonstrated [16] that there is a large region of the smoothing width where the 
shell correction is practically constant (plateau condition). Therefore the new curvature 
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correction methods have to be checked in the case of harmonic oscillator potential to 
see if the plateau is remained or not. 

The parameter of the harmonic oscillator potential {hw =6 MeV), the considered 
particle number {Nt-=70) and the order of the curvature correction (m = 3) were chosen 
to be the same as in p^. In this work the PCC (Gaussian type kernel) and the 
BPCC (Cosh-type kernel) methods were used. The shell correction as a function of 
the smoothing width is displayed in Figure [1] using Gaussian starting kernel. Since the 
shell spacings are hw the shell corrections are calculated only for 7 > /lo; otherwise the 
effect of shells can be noticed in the smooth SPLD [22] . 

The curve signed by PCC m = 3 in Figure [1] agrees with the result of [IE]. In 
order to see the dependence on the order of the curvature correction the calculation is 
repeated for m = 5 and it is displayed also in Figure 1. On the considered 7 range the 
plateau region of the MWCC method is roughly two times larger than the plateau of the 
PCC method. Furthermore the increase of the order of the curvature correction from 3 
to 5 does not modify the result of the MWCC method but there is large change in the 
case of the PCC method. 

In the case of Cosh-type starting kernel the shell correction is shown in Figure 
121 Using this kernel the quality of the DCC, MWCC and BPCC curvature correction 
methods can be judged. In this case also the MWCC method is proved to the best, it 
gives the largest plateau region. The results of the BPCC and MWCC methods differ 
only at large 7 values. In this large 7 region only the MWCC method keeps the plateau 
property. The importance of the large plateau region derives from the fact that in global 
mass calculation for a given nucleus only one 7 value is considered. It would be desirable 
to have this fixed 7 value in the plateau region. 

It remains to study the properties of the curvature correction methods in the case of 
realistic potential. We restrict ourselves for spherically symmetric mean field. We will 
consider Saxon- Woods potential with spin-orbit term. The geometrical and strength 
parameters of the Saxon- Woods potential are taken from the parametrization given in 
[17] . There is a great difference between harmonic oscillator and realistic potentials. 
The first one has only bound states whereas the Saxon- Woods potential supports both 
bound and scattering (continuum) states. 

It is known that the presence of continuum states and the finiteness of the number 
of bound states cause difficulties in the shell correction method [HI [18]. The formula ([1]) 
for the quantum SPLD contains only the contribution of bound states. The exact form 
of the continuum SPLD is known [18] and it was used in [T8| IT9] . The continuum SPLD 
is expressed by the help of the scattering phase shift. In shell correction calculations the 
bound states usually are obtained by diagonalizing the Hamiltonian on a given square 
integrable basis. On such a basis however it is difficult to calculate the continuum 
SPLD. An elegant way to solve this problem was developed in [20] and was used in 
shell correction calculation |21]. In this method two diagonalizations have to be carried 
out. The first diagonalization is the standard one (we take the Hamiltonian as it is). 
In the second diagonalization the interaction is switched off and only the kinetic energy 
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is diagonalized using the same basis as in the first diagonahzation. The results of the 
second diagonahzation are denoted by e.-^"* and the corresponding degeneracies by df'\ 
The smooth SPLD which contains the effect of scattering states reads 

*) = ^E(*^(^)-4°'a-(^)), (34) 

where the summation over i extends to all states gained from the diagonalizations. 

The value of the parameter c is determined in the same way as in the harmonic 
oscillator case. Since now we have only finite number of bound states the upper bound 
of the averaging window has a natural value X^S'i ^i- It corresponds to a nucleus where 
all bound states are occupied. The lower bound of the averaging window was N^- — 50. 
We have checked that the result is hardly depends on this value. 

The neutron shell correction for ^^^Gd is displayed in Figures 3 and 4. In the case 
of Gaussian starting kernel (Figure 3) the order of the curvature correction was three 
and five. Results using Cosh-type starting kernel are given in Figure 4 here the order of 
the curvature correction was seven. Prom Figures 3 and 4 the same conclusion can be 
drawn. In the case of the method MWCC a very large plateau region is present (the shell 
correction is practically independent from the value of 7). The other shell correction 
methods have no plateau, they satisfy only the local plateau condition (the derivative 
of the shell correction with respect to 7 is zero at some 7 point). The plateau property 
can be made quantitative if the total variation of the shell correction is determined in a 
given 7 region. Considering Figure 3 the total variation are 0.3 MeV and 3 MeV for the 
MWCC and the PCC methods. In the case of the harmonic oscillator the total variation 
of the shell correction is only 0.01 MeV for the MWCC method. The plateau property 
in the realistic potential case is not as good as in the harmonic oscillator case. 

The dependence on the order of the curvature corrections is the same as in the 
harmonic oscillator case. The results of the MWCC method is practically independent 
from the parameter m. In the case of the methods PCC, DCC and BPCC the shell 
correction is sensitive to m. 

5. Summary 

In the Strutinsky's shell correction method the smooth single particle level density 
is calculated by a convolution. The aim of the curvature correction is to satisfy a 
self consistency condition, the smoothing (convolution) must leave a smooth function 
unchanged. We are aware only two curvature correction methods but in global mass 
calculations only one of them is used. We call the standard method polynomial curvature 
correction since the curvature correction is achieved by multiplying the kernel of the 
convolution by a polynomial. Two new curvature correction methods are introduced. 
One of the new methods (we call it many width curvature correction (MWCC)) is proved 
to be the best among the investigated curvature correction methods. 
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For potentials without scattering states (harmonic oscillator and similar potentials) 
the shell correction is practically constant (plateau condition) on a very large region over 
the smoothing width in each curvature correction methods. However the length of the 
plateau in the MWCC method is much larger than in the other methods. 

For potentials supporting scattering states (realistic potentials) the advantage of 
the MWCC method is more pronounced. In the case of each methods except the MWCC 
method only the local plateau condition is satisfied. Due to the lack of convergence in 
the order of the curvature correction it is difficult to define a unique value for the shell 
correction. If the calculation is carried out by the MWCC method a large plateau region 
exists even in the case of realistic potentials and the shell correction is quite stable with 
respect to the increase of the order of the curvature correction. It seems to us that 
the MWCC method is better suited for shell correction calculation than the standard 
procedures. 
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Figure 1. Shell correction for the harmonic oscillator potential {Hui = 6 MeV, Nr~70) 
as a function of the smoothing width using different curvature correction methods. 
The order of the curvature correction is three and five. The kernels are derived from 
Gaussian starting kernel. In this case the methods PCC, BPCC and DCC give the 
same kernel for the smoothing. 
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Figure 2. Same as Figure 1 but the kernels are derived from Cosh-type starting 
kernel. The order of the curvature correction is three. 



Curvature Correction in the Strutinsky's Method 14 



> 



2 

-2 
-4 
-6 
-8 
-10 
-12 
-14 



"1 ' 

PCC m=3 
PCC m=5 
MWCC m=3 
MWCC m=5 




20 
[MeV] 



Figure 3. Neutron shell correction for ^™Gd as a function of the smoothing width 
using different curvature correction methods. The single particle mean field was a 

Saxon- Woods potential with spin-orbit term. The order of the curvature correction is 
three and five. The kernels are derived from Gaussian starting kernel. In this case the 
methods PCC, BPCC and DCC give the same kernel for the smoothing. 
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Figure 4. Same as Figure 3 but the kernels are derived from Cosh-type starting 
kernel. The order of the curvature correction is seven. 



